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Abstract. 

Statistical analysis of Dst behaviour during recovery phase of magnetic storms 
induced by different types of interplanetary drivers is made on the basis of 
OMNI data in period 1976-2000. We study storms induced by ICMEs (in- 
cluding magnetic clouds (MC) and Ejecta) and both types of compressed re- 
gions: corotating interaction regions (CIR) and Sheaths. The shortest, mod- 
erate and longest durations of recovery phase are observed in ICME-, CIR- 
, and Sheath-induced storms, respectively. Recovery phases of strong (Dst min < 
— 100 nT) magnetic storms are well approximated by hyperbolic functions 
Dst(t) = a/(l + t/rh) with constant times for all types of drivers 
while for moderate (—100 < Dst min < —50 nT) storms Dst profile 
can not be approximated by hyperbolic function with constant because 
hyperbolic time increases with increasing time of recovery phase. Rela- 
tion between duration and value Dst min for storms induced by ICME and 
Sheath has 2 parts: Dst m i n and duration correlate at small durations while 
they anticorrelate at large durations. 
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1. Introduction 

Dynamics of a magnetic storm is result of competing processes: excitation and relax- 
ation [Burton et al, 1975]. Generation of magnetic storm is connected with southward 
(Bz < 0) component of interplanetary magnetic field (IMF) (or interplanetary electric 
field Ey = Vx x Bz, where V is solar wind velocity) [Dungey, 1961; Fairfield and Cahill, 
1966; Rostoker and Falthammar, 1967; Russell et al, 1974; Burton et al, 1975; Akasofu, 
1981]. Relaxation of magnetic storm is connected with ring current decay [Daglis et al, 
1999]. The relation between excitation and relaxation changes during storm development: 
excitation process prevails over relaxation process at the main phase while the relaxation 
prevails at the recovery phase. 

As well known, large and long component Bz < in the interplanetary space is ob- 
served only in disturbed types of solar wind. Therefore interplanetary drivers of magnetic 
storms are the following disturbed types of solar wind: (1) CIR (Corotating Interaction 
Regions) formed in compress region between slow and high-speed streams of solar wind, 
(2) ICME (interplanetary coronal mass ejections including magnetic clouds and less pow- 
erful disturbance - Ejecta), and (3) Sheath formed in compress region between fast ICME 
and slower stream of solar wind (see reviews and recent papers, for instance, by Tsurutani 
and Gonzalez [1997]; Gonzalez et al. [1999]; Yermolaev and Yermolaev [2006]; Zhang et 
al. [2007]; Turner et al. [2009]; Yermolaev et al. [2010b, 2011]; Nikolaeva et al. [2011]; 
Gonzalez et al. [2011]; Guo et al. [2011] and references therein). One of recent important 
experimental results is evidence that features of magnetic storm depend on type of the 
interplanetary driver [Borovsky and Denton, 2006; Denton et al, 2006; Huttunen et al, 
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2006; Pulkkinen et al, 2007a; Plotnikov and Barkova, 2007; Longden et al, 2008; Turner 
et al, 2009; Yermolaev et al, 2010b; Guo et al, 2011] These facts indicate that mecha- 
nisms of magnetic storm generation (or modes of these mechanisms) can differ depending 
on the driver. 

It is usually considered that decay of the ring current depends exclusively on internal 
processes (charge exchange, Coulomb collision, wave-particle interaction, and drift loss) 
and current interplanetary conditions (see, for instance, Daglis et al. [1999]; Keika et al. 
[2006]; Xu and Du [2010]. and references therein) In this paper we study time variation of 
Dst index during recovery phase of magnetic storms induced by different interplanetary 
drivers, i.e. dependence of Dst profile on conditions which were in interplanetary space 
during the main phase. 

In accordance with formula by Burton et al. [1975] the Dst index should grow with 
exponential law at recovery phase. Many papers showed that for Dst profile approxima- 
tion it is necessary to use exponential time (r e ) which changes with time and depends 
on minimum Dst index [OBrien and McPherron, 2000; Feldstein et al, 2000; Monreal 
MacMahon and Llop-Romero, 2008; Xu and Du, 2010]. Aguado et al. [2010] indicated 
that Dst profile can be approximate by hyperbolic function Dst(t) = Dst /(1 + t/r h ) 
with constant hyperbolic time r h . In our analysis we use hyperbolic approximation as 
simpler one because our aim is to search for difference of Dst profile for storms generated 
by various interplanetary drivers and type of approximation function is not important for 
this aim. 
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2. Methods 

In order to study the problem we use OMNI data of solar wind and interplane- 
tary magnetic field parameters (see http://omniweb.gsfc. nasa.gov [King and Papi- 
tashvili, 2004]) and data on Dst index (see http://wdc.kugi.kyoto-u.ac.jp/index.html), 
as well as our catalog of large-scale interplanetary events for period of 1976-2000 (see 
ftp://www.iki.rssi.ru/pub/omni [Yermolaev et ai, 2009]). We study moderate and strong 
magnetic storms with —100 < Dst min < —50 nT and Dst min < —100 nT, respectively. 
The technique of determination of connection between magnetic storms and their inter- 
planetary drivers consists in the following. If the minimum of Dst index lies in an interval 
of a type of solar wind streams or is observed within 1-2 hours after it we believe that the 
given storm has been generated by the given type of streams [Yermolaev et ai, 2010a]. 
Our analysis for period of 1976-2000 showed that 145 magnetic storms (31.2% from total 
number of identified storms) have been generated by CIR, 96 (20.7%) storms by Sheath, 
62 (16.0%) by MC and 161 (34.7%) by Ejecta. The sources of other 334 magnetic storms 
(i.e., 42% of 798 storms observed during this time interval) are indeterminate and these 
storms are indicated as IND [Yermolaev et ai, 2010a, 2012]. About 20% of storms were 
multistep ones during recovery phase and these storms were excluded from analysis. 

It is should be noted that average duration of main phase of magnetic storm is about 
7 ± 4 hours and it is less than average durations of interplanetary drivers: 24 ± 11 h for 
MC, 29 ± 5 h for Ejecta, 16 ± 3 h for Sheath before Ejecta, 9 ± 5 h for Sheath before MC, 
and 20 ± 4 h for CIR [Yermolaev et ai, 2007, 2010a]. Therefore the main phase for the 
majority of magnetic storms may be completely defined by one type of solar wind streams 
and it is possible to consider that a storm is connected with one interplanetary driver. In 
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contrast with main phase, the recovery phase can last up to 3 days, i.e. it is longer than 
duration of interplanetary drivers. Therefore only at initial (short in comparison with 
full phase) stage of recovery phase the interplanetary conditions may correspond to that 
driver which operated at the main phase, and at the subsequent stage of recovery phase 
type of solar wind streams is replaced by another types: Fast follows CIR, MC or Ejecta 
follow Sheath, Fast or Slow streams follows MC and Ejecta. Change of type of solar wind 
during recovery phase is not analyzed in the given work, and we classify storms only on 
their drivers operating at the main phase of storm. 

One of the difficulties of recovery phase analysis is definition of time when the phase 
comes to its end as Dst index quickly grows at the first stage, and then the growth 
slows down at approach to the initial undisturbed level, and there is large (up to 20%) 
Dst variation which is not connected with general process of storm evolution. We use 
times when Dst index achieves levels of 1/2 and 1/3 from minimum Dst index as criteria 
of time of recovery phase termination, and analyze two durations Ati/ 2 and Ati/ 3 , i.e. 
time intervals from Dst mim up to (\j1)Dst mim (At 1/2 = t((l/2)Dst mim ) - t(Dst mim )) 
and (1/3) Dstmim (Aii/3 = t((l/3)Dst m i m ) — t(Dst mim )), respectively. Comparison of two 
data sets corresponding to Ati/ 2 and Ati/ 3 allows us to make conclusions about dynamics 
of Dst index during storm recovery phase. 

3. Results 

Histograms of recovery phases durations At±/ 2 (black line) and Ati/ 3 (gray line) for the 
storms generated by various drivers are shown in Figure 1. All distributions Ati/ 2 have 
peaks in range of 6-12 hours, some events in range of 0-6 hours and long tails in range 
of large times. In contrast with At 1 / 2 , character of distribution Ati/ 3 depends on type of 
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drivers. For Ejecta and MC (and their sum) distributions have a maximum in range of 
12-18 hours, approximately half of maximum in range of 6-12 hours and sharply falling 
down tail: ~ 1/10 of maximum for Ejecta and ~ 1/3 for MC at 42 hours. Histogram 
for CIR has a maximum in range of 24-30 hours and wide distribution. For Sheath 
there are 2 maxima at 12-24 hours (this maximum is similar to maximum Ejecta and 
MC ) and 36-42 hours and wide distribution. IND-storms (with not identified sources) 
have the distribution similar to the sum of mentioned above distributions. Mean values 
and standard deviations of durations of recovery phases for both Atx/ 2 and Ati/ 3 and 
for all types of drivers (see Table 1) quantitatively support features of Figure 1. Thus, 
both figure and table show that distributions differ for ICME and compressed types of 
drivers. At transition from Ati/ 2 to Atx/ 3 distributions for Ejecta and MC (and their 
sum) are simply displaced towards longer durations with distribution form preservation 
(sharp maximum at small durations and monotonously falling down tail in the range of 
large durations), while CIR- and Sheath-induced storms qualitatively change the shape 
of distributions (larger time of maxima and wider disributions). 

Figure 2 shows Dst profiles during recovery phase of strong (Dst m i n <-100 nT) mag- 
netic storms induced by various drivers, and their approximations by hyperbolic functions 
(approximation coefficients are presented in Table 2). Data for every type of drivers were 
obtained by method of superposed epoch analysis with zero (epoch) time equal to time of 
Dstmim. All profiles are well approximated in region of 0-48 hours, the largest time is 
observed for Sheath-induced storms and the smallest for ICME-induced storms. Unlike 
strong storms, moderate storms (see Figure 3) for all drivers have a characteristic break in 
range of 10-15 hours, and they cannot be described by hyperbolic function with constant 
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Th (or exponential function with constant r e ). For example, we show in Figure 3a, 3b and 
3c results of approximation in ranges of 0-11, 0-18 and 0-48 hours, respectively, and we 
present approximation coefficients in the Table 2. For all approximation regions time 
is the largest for Sheath-induced storms and the least for ICME-induced storms. 

In order to investigate how duration of storm recovery phase depends on value Dst min , 
we calculate average \Dst m i n \ values for each type of drivers in 6-h bins, and results 
are shown in Figure 4. Panels 4a and 4b show results for durations Aii/ 2 and Ati/ 3 , 
respectively. For all types of storms on the average the size of storms grows in a range 
from 4 till 18 hours for Ati/ 2 durations and in a range of 6-30 hours for Ati/ 3 , and falls 
with increasing duration. For At 1/3 the size of storms is practically a constant for duration 
longer 48 hour, i.e. it does not depend on duration. Near maxima \Dst m i n \ stronger storms 
were induced by Sheath and MC, Ejecta-induced storms have moderate peak and CIR- 
induced storms have weak maximum. It should be noted that maxima \Dst min \ on both 
panels are statistically significant (see Figure 1). Thus, there are 2 branches of Dst min 
versus duration dependence: increasing and decreasing ones. 

4. Discussion and conclusions 

On the basis of OMNI data of plasma and IMF parameters of solar wind during 1976- 
2000 we classified various types of solar wind streams and found interplanetary drivers for 
572 magnetic storms. These data allowed us to compare temporal evolution of Dst index 
during recovery phase of magnetic storms induced by CIR, Sheath and ICME (including 
MC and Ejecta). Our study allowed to obtain following results. 

I. As distributions of duration of recovery phases show (see Figure 1 and Table 1), 
initial parts of recovery (distributions of Ati/ 2 ) are close to each other for all types of 
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drivers, and at further part of recovery (Ati/ 3 ) distributions differ and histograms are 
wider (i.e. on the average recovery is slower) for compression regions Sheath and CIR 
than for both types of ICME. 

2. Dst index during recovery phase is well approximated by hyperbolic functions 
Dst(t) = a/(l + t/r h ) only for strong (Dst min < —100 nT) storms. This result is in 
good agreement with conclusions by Aguado et al. [2010]. Approximation shows that the 
hyperbolic index is the greatest for Sheath- (the slowest recovery), intermediate for 
CIR- and the least for ICME-induced storms (Figure 2 and Table 2). 

3. In contrast to strong storms, Dst index during recovery phase is bad approximated by 
hyperbolic functions for moderate storms with —100 < Dst m i n < —50 nT and hyperbolic 
index depends on time of recovery phase. Approximation in regions of 0-11, 0-18, and 
0-48 hours (see Figure 3) shows that dependence of hyperbolic index on type of drivers 
is similar to dependence for strong storms. 

4. For storms generated by ICME and Sheath, there are 2 classes of storms (see Figure 
4): duration of recovery phase and Dst min correlate for short storm duration and they 
have inverse correlation for long storm duration. For storms generated by CIR, duration 
of recovery phase of storm does not depend on Dst min . 

Thus obtained results allow us to conclude that recovery phase Dst variations depend 
on type of interplanetary drivers inducing magnetic storms and magnetosphere remember 
type of driver during recovery phase. 
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Table 1. Average values and standard deviations of recovery phase durations. 
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Ati/2 = 



t(l/2Dst min ) - t(Dst min ) 

b At 1/3 = t{l/3Dst min ) - t(Dst min ) 
Table 2. Approximation of Dst profiles by hyperbolic function Dst(t) — a/(l + t/rh) 
for moderate and strong magnetic storms induced by various drivers 



Types of driver 


-100 < Dst min < -50 


Dst min < -100 


(0 - ll)h* 
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-68.4 14.97 
-66.4 15.43 
-70.1 13.39 
-71.5 16.26 
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-66.3 18.62 
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-69.6 14.05 
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-59.2 35.46 
-59.2 29.94 
-58.8 32.36 
-60.0 42.74 
-59.0 34.72 
-59.5 48.31 


-141.5 17.95 
-144.1 19.69 
-135.1 15.72 
-154.3 16.95 
-142.8 16.23 
-148.8 22.37 



* Time interval of approximation. 
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Figure 1. Distributions of magnetic storm recovery durations At\/2 (black line) and Ati/3 



(gray line) for different types of interplanetary drivers 
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Figure 2. Temporal evolution of recovery phases of strong (Dst m i n < —100 nT) storms 
induced by various types of drivers (points) and hyperbolic approximations in time intervals of 
0-48 hours. 
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Figure 3. Temporal evolution of recovery phases of moderate (—100 < Dst min < —50 
nT) storms induced by various types of drivers (points) and hyperbolic approximations in time 
intervals of (a) 0-11 hours, (b) 0-18 hours, and (c) 0-48 hours. 
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induced by various types of drivers: (a) for Aii/ 2 and (b) for Ati/ 3 
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